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Axions are very light, very weakly interacting particles postulated more than 30 years ago in the 
context of the Standard Model of particle physics. Their existence could explain the missing dark 
matter of the universe. However, despite intensive searches, they have yet to be detected. In this 
work, we show that magnetic fluctuations of topological insulators couple to the electromagnetic 
fields exactly like the axions, and propose several experiments to detect this dynamical axion field. 
In particular, we show that the axion coupling enables a nonlinear modulation of the electromagnetic 
field, leading to attenuated total reflection. We propose a novel optical modulators device based on 
this principle. 



PACS numbers: 75.30.-m, 78.20.-e, 78.20.Ls, 03.65.Vf 

The electromagnetic response of a three dimensional 
insulators is described by the Maxwells action 5*0 = 
^ / d? xdi^eEi"^ — ^B^), with material-dependent dielec- 
tric constant e and magnetic permeability /x, where 
E and B are the electromagnetic fields inside the in- 
sulator. However, generally, it is possible to include 
another quadratic term in the effective action Sg — 
Itt'Tk J d^xdtEi ■ B, where a = /he is the fine struc- 
ture constant, and is a parameter describing the in- 
sulator in question. In the field theory literature this 
effective action is known as the axion electrodynamics^, 
where plays the role of the axion field. Under the pe- 
riodic boundary condition the partition function and all 
physical quantities are invariant if 9 is shifted by integer 
multiples of 27r. Therefore all time reversal invariant in- 
sulators fall into two distinct classes described by either 
9 = {) OY 9 = Ti^ . Topological insulators are defined by 
9 = TT and can only be connected continuously by time re- 
versal breaking perturbations to trivial insulators defined 
by 9 = 0. The form of the effective action implies that an 
electric field can induce a magnetic polarization, whereas 
a magnetic field can induce an electric polarization. This 
effect is known as the topological magneto-electric effect 
(TME) and 9 has the meaning of the magneto-electric 
polarization P3 = 9/2Tr. Physically the parameter 9 de- 
pends on the band structure of the insulator and has a 
microscopic expression of the momentum space Chern- 
Simons form^ 
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where A^^Ck) — — i(Q!k|^-|/3k) is the momentum space 
non-abelian gauge field, with indices a,P referring to 
the occupied bands. The 9 parameter has been calcu- 
lated explicitly for several basic models of topological 
insulators^''^. In a topological insulator the axion field 
gives rise to novel physical effects such as the image 
monopole and anyonic statistics'*. This field, however, 
is static in a time-reversal invariant topological insula- 
tor. In this work, we consider the anti-ferromagnetic long 
range order in a topological insulator, which breaks time- 



reversal symmetry spontaneously, so that 9 becomes a 
dynamical axion field taking continuous values from to 
2tt. In the following we will refer to such an antiferromag- 
netic insulator as a "topological magnetic insulator" . We 
propose a minimal model in which the antiferromagnetic 
order break the time reversal symmetry spontaneously 
and the magnetic fiuctuations couple linearly to the axion 
field, thus realizing the dynamic axion field in condensed 
matter systems. Compared to its high energy version, 
the axion proposed here has the advantage that it can be 
observed in controlled experimental settings^. With an 
externally applied magnetic field, the axion field couples 
linearly to light, resulting in the axionic polariton. By 
measuring the attenuated total refiection, the gap in the 
axionic polariton dispersion can be observed. An attrac- 
tive feature is that the axionic polariton gap is tunable 
by changing the external electric or magnetic fields. The 
control of the light transmission through the material 
enables a novel type of optical modulator. We also pro- 
pose another experiment to detect the dynamic axion by 
microcantilever torque magnetometry, where the double 
frequency response of the cantilever is a unique signature 
of the dynamic axion field. 

We propose several materials that may realize the 
topological magnetic insulator with dynamic axion field. 
One possibility is the topological insulator Bi2Te3, 
Bi2Se3, Sb2Te3 doped with 3c? transition metal elements 
such as Fe^ii. Another possible class of material is the 
5d transition metal compound AxBOy with B and A 
standing for some 5d transition metal and some alkali 
metal, respectively. Electrons in 5d-orbital can have both 
strong spin orbital coupling and strong interaction, which 
is ideal for the realization of the topological magnetic 
insulator—. We show that such a compound with the 
corundum structure may have a topological magnetic in- 
sulator phase if the states closed to fermi level are formed 
by t2g orbitals with total angular momentum Joff = 1/3^. 
The detail of this proposal is beyond the goal of the 
present work, and will be presented in a separate paper.— 
We also noticed two very recent works on 5d transition 
metal compounds with pyrochlore structure, which may 
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also realize the topological magnetic insulator phaseJ^^ 
Effective model for the 3D topological insulator Al- 
though all the physical effects discussed in this paper are 
generic for any system supporting axionic excitation and 
do not rely on a specific model, we would like to start 
from a simple model for concreteness. We adopt the effec- 
tive model proposed by Zhang et al in Refi^ to describe 
topological insulators Bi2Te3, Bi2Se3 and Sb2Te3. The 
low energy bands of these materials consist of a bond- 
ing and an anti-bonding state of Pz orbitals, labeled by 
|P2-,t (i)) and |P1+,T (i)), respectively. The generic 
form of the effective Hamiltonian describing these four 
bands is obtained up to quadratic order of momentum k 
in Reffi^. Since a lattice regularization is necessary for 
computing axion field 6, in the present paper we start 
from a lattice version of this model, with the Hamilto- 



Hoik) = eo(k)l4x4 + ^da(k)r'^ 



(2) 



di,2,...,5(k) = {A2Smkx:,A2sinky,Aismkz,M{'k),0) 

where eo(k) = C + 

2Di + AD2 - 2Di cos kz - 2D2 (cos k^ + cos ky), M (k) = 
M - 2Bi - AB2 + 2Bi cos kz + 2B2 (cos k^ + cos ky), and 
the Dirac T matrices have the representation r*^^'^''^''*'^) — 
(ax 'S) Sx,(7x Sy,ay ig) I2x2,crz <S) I2x2,crx ® Sz) in the 
basis of {\Plt, T), \Plt, i), \P2:A), \P27, i)). 

V,T breaking terms. The above Hamiltonian pre- 
serves both time reversal symmetry T and parity V. 
In the effective model ([2]), the time- reversal and 
spatial inversion transformation are defined as 
T = i{l2y.2 ® Sy) ■ IC (with K the complex conju- 
gation operator) and V = cr^ ® 12x2, respectively. The 
axion field 9 can be calculated by formula ^ which 
gives or TT depending on the value of parameters^, 
as expected from time reversal symmetry. Now we 
consider a perturbation to the Hamiltonian ([2]) which 
can lead to deviation of 9 from or tt. Since 9 is 
odd under time reversal and parity operation, only 
time reversal and parity breaking perturbations can 
induce a change of 9. By simple algebra one can show 
that to the leading order the T ,V breaking perturba- 
tion must have the form SHik) = X]i=i 2 3 5 
Thus the perturbed Hamiltonian can still be 
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written as i/(k) c^a(k)r" with da(k) = 

a=l 

{A2 sin kx + mi, A2 sin ky + 7712, Ai sin k^ 4- m^, Al(k), TO5). 
For this model, Eq. ^ can be reduced to an explicit 
expression for 9: 
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d5 {dk^ di ) {dky d2 ){dk,d3). 

(3) 

Although there are four independent parameters TOi.2,3,5 
in T, V breaking term, only leads to a correction to 
9 to the linear order. Thus in the following we will take 
'1^1,2,3 = without the loss of generality, and leave only 
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FIG. 1: [Color Online] Crystal structure of Bi{Fe)2Sea with 
three primitive lattice vectors denoted as ti,2,3. A quintuple 
layer with Sel - Bi{Fe) 1-Se2- Bi{Fe) 1' - Sel' is indicated 
in the orange box. The spin ordering configuration giving 
rises to the Fs mass is indicated by the black arrow, which is 
antiferromagnetic along z direction and ferromagnetic within 
xy plane. 



the TOsF^ term. To see the physical meaning of the msF^ 
term we change the basis to {\A, f), \A, [), \B, t), \B, [)) 
with lAT (i)) = ;^(l^l+,T(i)) + l^2-,t(i))) and 

\BA (i)> = j^mtA a))-\P2-,Hl)))- We see 
that V\AA (i)> = |i?,T (i)). Physically, |A,T (i)) and 
\Bi T (i)) stand for states of the two lattice sites in each 
unit cell, e.g. the two Bi atom sites in Bi2Te3, which 
are shifted away from the inversion center along z di- 
rection. By transforming F^ to the new basis, we see 
that it represents a staggered Zeeman field that points in 
+z{—z) direction on the A{B) sublattice. This staggered 
Zeeman field can be generated by an antiferromagnetic 
order, where electron spins point along opposite z direc- 
tions on two sublattices, as shown in Fig. 1. Without 
the F^ term, topological insulators have protected sur- 
face states consisting of odd number of massless Dirac 
cones. With broken time reversal symmetry, a direct 
consequence of the m^T^ term is that it opens a gap in 
the surface state spectrum. The surface state gap is equal 
to m5, independent on the orientation of the surface. 

Dynamical axion in the anti-ferromagnetic phase We 
have seen above that the axion field 9 can deviate from 
or TT if the electrons are coupled to the antiferromag- 
netic order parameter. Such a staggered field can be 
induced if an antiferromagnetic long range order is es- 
tablished in the system. Physically, an antiferromagnetic 
ordered phase can be obtained naturally if the screened 
Coulomb interaction of electrons become strong. For ex- 
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ample, in the Bi2Te3 family doped with magnetic impu- 
rities such as Fe, the substitution of p-electrons of Bi or 
Te by d-electrons of the magnetic impurities effectively 
enhanced the on-site repulsion of electrons. In the basis 
(|A,t),|A, i),|B,t),|B,i)) the Hamiltonian with inter- 
action can be written as 

H = Ho + {niA^riiAi + niB^riiBi) + V'^niAniB- 

i i 

, (4) 

in which the first term is the kinetic energy given by 
Eq. ^ and the rest two terms represent the on-site 
repulsion U and the inter-site repulsion V between A 
and B sites. Possible ordered phase resulting from the 
interaction include the ferromagnetic phase where the 
spin on two sublattices A and B point in the same di- 
rection, the antiferromagnctic phase where the spin on 
two sublattices point in the opposite direction, and the 
charge density wave (CDW). Correspondingly the order 
parameters are taken as the ferromagnetic order param- 
eter M+ — \ {{^iA) + (Sis)), the antiferromagnctic or- 
der parameter = \ {{^iA) - (Sis)) and the CDW 
order parameter p = ^ {{uiA) — (nis)). It is assumed 
that translational symmetry is preserved and all the or- 
der parameters are uniform in space. In the mean field 
approximation, we find that for a wide range of value for 
band structure parameters M, Ai^2 and B12, the sys- 
tem develops antiferromagnctic order pointing in the z 
direction 

of V, which thus leads to TO5 

Axion electrodynamics. In the mean-field approxima- 
tion, the antiferromagnctic phase has a static axion field 
9. However, the antiferromagnctic phase also has am- 
plitude and spin wave excitations, which can induce 
fluctuations of the axion field. The fluctuation of the 
Neel vector M~ can be generally written as M~ = 
(Mp" + 6M^{x, t)) z + SM^{x, t)x + 6My{x, t)y. To the 
linear order, it can be shown from symmetry analysis 
that the fluctuation of axion field only depends on 5Mz, 
since 6* is a pseudo-scalar. In other words, we have 
(50(x, t) = 6m5{-K,t)/g — — '^U5M.,{yi,t)/g where the co- 
efflcient g can be determined from Eq. The disper- 
sion of the amplitude mode 5Mz (x, t) can be obtained in 
standard RPA approximation, leading to a massive ax- 
ion field (5^?(x, t). Considering the coupling term 6'E • B 
of axion with electromagnetic field, the effective action 
describing the axion-photon coupled system is given by 



A/q z if the effect of U dominates that 
- 1 [/ Mr and axion field 
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d^xdticE^ - -B^) 

IT J fi 

+ ^ J d^^dt {9o + (561) E • B 



trie field and the magnetic field respectively, e and fj, 
are the dielectric constant and magnetic permeability 
respectively. The second term describes the topologi- 
cal coupling between the axion and the electromagnetic 
field, with a = /he the fine-structure constant. The 
third term describes the dynamics of the massive axion. 
Within the model we have adopted, the parameters J 
and m are given by 
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(27r)3 4|d|3 



(6) 



where \d\ — yJ2a=idad'^ and the repeated index indi- 
cates summation with i = 1, 2, 3, 4. 

The axionic polariton. The dynamic axion field 9 
couples nonlinearly to the external electromagnetic field 
combination E • B. When there is an externally applied 
static and uniform magnetic field Bq parallel to the elec- 
tric field E of the photon, 9 will couple linearly to E^^. In 
condensed matter systems, when a collective mode is cou- 
pled linearly to photons, hybridized propagating modes 
called polariton emergei^. The polaritons can be coupled 
modes of optical phonon and light through the electric 
dipole interaction, or coupled modes of magnon and light 
through the magnetic dipole interaction. Here we pro- 
pose a novel type of polariton — axionic polariton which 
is the coupled mode of light and the axionic mode of an 
antiferromagnet. The dispersion of the axionic polariton 
can be obtained from the effective action ([5]) which leads 
to the following linearized equation of motion- '-^^ 
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,T o aBn d 
E - c'^V^E + - ° 
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^59~v''y^59 
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Sn^g^J 



= 



E = 0, (7) 



where c' is the speed of light in the media and e is the 
dielectric constant. Compared to photon, the dispersion 
of axion can be neglected, in which case the axionic po- 
laritons have the dispersion 



± ]^\l{c''^k'^ +w? + &2)2 _ 4c'2fc2m2, (8) 

with 52 — a^Bo' I'^TC'cg^ J . As shown in Fig. 2 a, this 
dispersion spectrum consists of two branches separated 
by a gap between m and \J m? -\- IP- . The quantity h mea- 
sures the coupling strength between the axion field and 
the electric field and is proportional to the external mag- 
netic field Bq. Upon turing on Bq, the axionic mode at 



d^xdt \{dtS9Y — {vidi59Y — m^50^] (5) k = changes its frequency from m to ^/ra^ + 62, due 



where J, Vi and m are the stiffness, velocity and mass 
of the spin wave mode 5Mz, E and B are the elec- 



to the linear mixing between the axion and the photon 
field. Physically, the axionic polariton is very similar to 
the transverse optical phonon polariton, since the axion 
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FIG. 2; [Color Online] Axionic polariton and ATR experi- 
ment, a. The dispersion of the axionic polariton. The grey 
area indicates the forbidden band between frequencies m and 
\Jm? + &2 (see text), within which light can not propagate 
in the sample. The red dashed line shows the bare photon 
dispersion uj — c'k. b. Setup for the attenuated total re- 
flection (ATR) experiment. Without external magnetic fleld, 
the incident light can transmit through the sample, c. When 
an external magnetic fleld is applied parallel to the electric 
fleld of light, the incident light will be totally reflected if its 
frequency lies within the forbidden band. 



also leads to an additional contribution to the charge po- 
larization due to the topological magneto-electric effect^, 
P = a'So/irO + eE. The optical phonon polariton has 
the same dispersion as Eq. ([5]), with the parameter b 
replaced by the lattice unscreened plasmon frequency 

Up = Y^ ^"^r^ ■ The key difference between axion and 
optical phonon is that the coupling between axion and 
electric field is determined by the external magnetic field 
Bo, which is thus tunable. 

The gap in the axionic polariton spectrum may be ex- 
perimentally observed using the attenuated total reflec- 
tion (ATR) method. The geometry is arranged such that 
the incident light is perpendicular to the surface of the 
sample and the static magnetic field is parallel to the 
electric field of light, as shown in Fig. 2 b and c. Since 
the light can only propagate through the media in the 
form of the axionic polariton, when the frequency of the 
incident light is within the gap of the axionic polariton 
spectrum, a significant increase of the reflectivity will 
be observed. To estimate the gap we adopt the tight- 
binding parameters obtained for 612863^'^ in Hamiltonian 
([2]). We take a typical exchange splitting for an antifer- 
romagnet 7715 = 15meV, and an estimated dielectric con- 
stant e — 100. With a magnetic fleld Bq = 2T, we obtain 
the axion mass m = 36.7meV and b = 2.2meV. The gap 
is approximately ^^^^ = 0.25meV, which can be observed 
experimentally. One unique signature of the axionic po- 
lariton is the dependence of the gap on Bq, which can be 
used to distinguish from usual magnetic polaritons. 
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FIG. 3: [Color Online] Cantilever Torque Magnetometry 
measurement of axion. A cantilever experiences a torque r 
when the magnetization Mt of the deposited topological mag- 
netic insulator induced by the AC electric fleld E(i) is placed 
in the applied magnetic fleld B. The u and 2a; frequency re- 
sponse comes from normal contribution and dynamical axion 
contribution, respectively. "iOSCAR" is so called interrupted 
oscillating cantilever-driven adiabatic reversal protocol which 
can be used to measure the frequency change of cantileveiii^. 



By changing the magnitude of Bq we can selectively de- 
termine the frequency band within which the light is to- 
tally reflected. This principle may flnd application as 
an amplitude optical modulator working at far-infrared 
frequency. 

Measuring the axion by microcantilever. As discussed 
in Refi^, a static 9 leads to a topological magneto-electric 



effect, in which a magnetization Mj = 



rSE is in- 



duced by electric fleld. When the dynamics of axion 
field 9 is considered, the behavior of the magnetization 
Mt is modified, which can be detected by microcan- 
tilever torque magnetometry (MTM).— li^i^S As shown 
in Fig. 3, a DC magnetic field B and an AC electric field 
E(i) = EiAc sin(tjt) are applied to the topological insu- 
lator attached to the tip of the cantilever, with an angle 
if between them. A magnetization Mt will be induced 
along the direction of the electric field, so that a torque 
T — Mt X B acting on the cantilever is generated. This 
magnetic force on the cantilever mimics a change in can- 
tilever stiffness, in turn, shifts the cantilever frequency by 
A/ = ^MtB^}Pi9^^ Here k is the cantilever stiffness, 

/o is the cantilever natural frequency, IcS is the effective 
length of cantilever. With the electric field E^c'Sin(a;t), 
the time dependent frequency change is, 



^fAc{t) — a\ sin(a;t) + a-i cos(2wi), 



(9) 



with 



fo 



47r2 2fc(U)2 



B9qEac sin ip, 



n2 r?2 



sin 2lp 



327r4 2fc(U)2 4g2 J(m2 - u"^) 



(10) 



It should be noticed that the response 02 with doubled 
frequency 2lo is a unique signature of the axion dynamics, 
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since all the response should have been linear and thus in 
frequency of lo if the axion field 6 is static. With typical 
cantilever parameters in present experiment o^^'^'^ /o = 
5.3kHZ, fc = 2.0 X lO-^N/m, U = SOOAim, B = lO^G, 
Exc = lO^'V/m, uj/2n — 1.5Hz and (p — 7r/4, we obtain 
the frequency shift ai — 120Hz and 02 = 4.04Hz, which 
are easily detectable in current experimental techniques. 

In conclusion, we have proposed the existence of dy- 
namical axions in topological magnetic insulators. We 
have proposed two experiments in which axions can be 
detected by its unique coupling to the electromagnetic 
field. The coupling between axion and photon leads to 



an axionic polariton which has a polariton gap tunable by 
magnetic field and thus may be used as a tunable optical 
modulator. 
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